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Abstract: We show that a version of the nearly Minimal Supersymmetric Standard Model 
(nMSSM), extended only in the singlet sector to include the additional superfields of right- 
handed neutrinos and very heavy Dirac particles conserving B — L, admits a viable scenario 
for Dirac leptogenesis and naturally small Dirac neutrino masses. The origin of the (B — L)- 
conserving high singlet neutrino scale and the desired supersymmetry breaking terms is 
associated with dynamical supersymmetry breaking in the hidden sector. 
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1. Introduction 

Leptogenesis, generating a sphaleron-induced baryon asymmetry from lepton asymmetry, 
can be realized not only with Majorana neutrinos |Q] but also with Dirac neutrinos conserv- 
ing B — L ||. Such neutrinos (as well as charged leptons) are first produced in the early 
universe with equal numbers of left-handed and right-handed states from the decays of very 
heavy fields. However, the Yukawa couplings of the former are small enough to hide, during 
the electroweak transition epoch, the right-handed lepton number from sphalerons. This 
is as opposed to the left-handed lepton number which gets converted to baryon number 
through sphaleronic transitions. 

In the above context, a realistic supersymmetric model for Dirac neutrino masses and 
leptogenesis was proposed in Ref. ||, and studied in detail in Refs. j|. In this model, a 
sub-eV neutrino mass scale gets generated by a combination of reasonably small Yukawa 
couplings and a high mass scale of heavy SU(2) doublet superfields. Furthermore, lep- 
togenesis can successfully arise from the decays of such heavy fields. However, there is a 
serious problem with this kind of a supersymmetric scenario. The same heavy superfields 
can also decay into right-handed sneutrinos which have masses of the order of hundreds 
of GeV to a bit more or less than a TeV, as induced by soft SUSY-breaking terms. Since 
their couplings yLHN, which are supersymmetric generalizations of Yukawas, are tiny 
with y ~ 1CP 13 guaranteeing the non-equilibration of the right-handed neutrinos, they do 
not decay quickly enough and persist as relics after the electroweak transition. In fact, 
their number density divided by the entropy density then is expected to be as large as 
10~ 3 . If they were stable, their relic density today would be far too much in excess of that 
expected of Dark Matter (DM) as estimated from cosmological observations. On the other 
hand, because of their tiny couplings, unstable right-handed sneutrinos can only decay 
very slowly, thus surviving till long after the decoupling of the Lightest Supersymmetric 
Particle (LSP) which takes place typically at a temperature of the order of 10 GeV. Their 
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late decays will then lead to a late non-thermal over-production of the LSPs, once again 
conflicting with the required DM relic density. Our aim is to present the conditions on 
supersymmetry breaking parameters in the neutrino sector avoiding this difficulty, and we 
present here a variant of the Murayama-Pierce model Q which achieves this goal. For 
non-supersymmetric models for Dirac leptogenesis, see Ref. ||. 

Our spectrum extends from that of the Minimal Supersymmetric Standard Model 
H only in the sector that is a singlet under the Standard Model gauge group. This 
extension is in the spirit of the nearly Minimal Supersymmetric Standard Model nMSSM 
Q. In the latter, one adds to the MSSM superpotential just one term, coupling a singlet 
chiral superfield S to the two Higgs superfields of the MSSM, in order to generate the 
//-term. In so far as direct effects on the low energy superpotential are concerned, we 
need to add to the nMSSM spectrum the right chiral neutrino superfields Ni, i standing 
for flavor, heavy Dirac pairs of singlet chiral superfields §k,^% carrying B — L charges 
±1, and a generic heavy singlet superfield X mediating supersymmetry breaking. In order 
to quicken the decays of relic right-handed sneutrinos, we need to have large values of 
certain supersymmetry breaking parameters whose contributions to the scalar potential 
get scaled by the neutrino mass and remain small without much affecting the lightest 
sparticle spectrum. Such large supersymmetry breaking parameters may arise consistently 
from a dynamical supersymmetry breakdown in the hidden sector. We show this explicitly 
by constructing a specific scheme in which the effects of such a breaking are fed down to 
the observable (neutrino) sector by the field X associated with U(l) m mediation ||, |9|. 

In Section II we describe our mechanism for Dirac neutrino leptogensis with a toy 
superpotential. In Section III we discuss the required strengths of the supersymmtery 
breaking parameters to avoid the late non-thermal overproduction of LSPs. The origin of 
such parameters from a dynamically broken supersymmetry in the hidden sector, fed down 
by U{\) gauge mediation, is outlined for a particular scheme in Section IV. Finally, Section 
V summarizes our conclusions. 



2. Mechanism for Dirac neutrino leptogenesis 

Consider first the toy superpotential 

WW = h ik LiH 2 $l + h' jk NjS^ k + Affc$ fe $£, (2.1) 

which will later constitute part of our full superpotential for the observable sector. Here 
h, h! are Yukawa coupling strengths and M k are heavy masses. We have used the standard 
notation for the superfields Li, Ni, H 2 , etc. Integrating out the heavy pairs, we have 

W eff = -JI 11 LiH 2 NjS = LiHiNjS, 2.2 
M k v 2 v s 

generating the Dirac neutrino mass matrix m"j with the vacuum expectation values (VEVs) 
V2 = {H®) an d v s = ($)■ Note that we need more than two pairs of heavy fields $?) in 
order to generate at least two non-vanishing neutrino mass eigenvalues. In this approach, 
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a tiny Dirac neutrino coupling y = m v jv2 ~ 1CP 13 is obtained from reasonable values of 
the Yukawa coupling strengths hh' ~ 10~ 6 and the heavy mass scale Mi ~ 10 10 GeV with 
vs ~ TeV. At this point we do not specify the origin of the VEV and the mass associated 
with the singlet superfield S. 

Let us assume two pairs of heavy superfields for simplicity, which gives us two non- 
vanishing neutrino mass eigenvalues. When the masses of the heavy neutrinos are hierar- 
chical, Mi <C M2, the decays of the lightest superfields <3?i, $f will produce the asymmetries 
e and — e in the final states LH 2 and NS, respectively. The decay rate of $1, is 

1 

8^ ^ 



r ^ = ^E[ 2 i^i 2 + i^ii 2 ] M i ( 2 - 3 ) 



and the asymmetry e, defined by eT D = r($i, $f — » LH 2 ) - r($i, $f — > LH 2 ), is 



(2.4) 



where 5 = M1/M2 @j- Here we assume that the heavy masses Mj are real and positive 
without loss of generality. For a successful leptogenesis, the lepton asymmetry, normalized 
by the entropy density s, Yl = (m — n l)/ s i 1S required to be 

Y L « e y$ ~io- 10 , (2.5) 

where Y$ is the out-of-equilibrium density of the heavy superfield pair <3?i,<I>f. This is 
usually expressed in terms of the ratio between the relativistic equilibrium number density 
and the entropy density j s ~ 10 -3 as well as the efficiency factor 77 < 1: Y$ = {n^ js)r\. 
The efficiency factor, which is controlled by the predominance of the decay over the inverse 
decay process, is determined by the ratio between the decay rate and the expansion param- 
eter of the universe H(T) at T = Mr. K = w "172 ^ Efc( 2 l /l fci| 2 + l /i fcil 2 ) where M P 
is the Planck mass and g* is the relativistic degree of freedom at T = M\ . The approximate 
functional form of the efficiency factor can be expressed as 7] ~ 1 /K (In K) 0,8 in the case of 
K ^> 1 [jn]]. From the relation (2.2), one typically gets K ~ (m l/ /10~ 3 eV)(TeV/f5) leading 
to K ~ 10 (and thus 77 ~ 0.1) for the solar neutrino mass scale m v ~ 0.01 eV and vs ~ 
TeV. A more specific choice of parameters for this will be shown below. To achieve the 
maximal efficiency tj ~ 1 corresponding to K ~ 1, we may take vs ~ 10 TeV or introduce 
one more pair of heavy fields associated with a smaller neutrino mass, m u ~ 10" 3 eV, as 
in the case of the usual seesaw mechanism [|10| . 

Let us remark that the asymmetry (2.4) is suppressed by the factor 5 <C 1 for hier- 



archical heavy neutrino masses. Nevertheless, a realistic value of e can be attained []|, 10] 
by choosing somewhat large magnitudes for h, h'. On the other hand, the CP asymmetry 
can be resonantly enhanced if S s» 1. This possibility with h ~ b! may well be motivated 
to generate the observed mild hierarchy in the light neutrino masses. In this case, the CP 
asymmetry is generalized to 

V 4 ^( 2 M 2 + K,l 2 )i-* 2 ' [ j 
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taking into account the contributions from the decays of two pairs $/,<5^ with I = 1,2. 
The above formula can be straightforwardly generalized to the case of three families of 
heavy superfields. For a detailed analysis of the parameter region satisfying the leptoge- 
nesis conditions, we refer the readers to Refs. [||, |10| ]. As an example of the parameters 
accommodating a successful leptogenesis, let us choose M, ~ 10 10 GeV and h ~ b! ~ 1CT 3 . 
This gives us K ~ 10 leading to rj ~ 0.1 and thus ~ 10" 4 M. Now taking 1 - £ ~ 0.1, 
we get the desired value of e ~ 10~ 6 . Note also that one can have almost degenerate heavy 
fields, 5 ~ 1, in which case the condition (2.5) can be met by significantly enhanced e and 
much smaller Y$. 

Independently of the details of the parameter space, one can conclude that successful 
leptogenesis via either (2.4) or (2.6) implies Y$ > 10~ 10 . This will lead in a supersymmetric 
scenario to the problem of unwanted relics which we propose to solve below. 



3. Conditions on supersymmetry breaking parameters 

A key requirement of Dirac leptogenesis is not to equilibrate the right-handed neutrinos 
with left-handed ones both above and at the electroweak phase transition when sphaleron 
interactions are active. This is needed to make the asymmetry in the left-handed neutrino 
turn into the baryon asymmetry via those interactions ||]. In other words, the scattering 
rates T$, induced by the effective operator (2.2), needs to be suppressed as compared 
with the expansion rate H of the universe above and at the electroweak phase transition: 
Ts < H for T > T c , where T c ~ 100 GeV. For the Dirac neutrino Yukawa operator 
W = (m u I '1J2) LH2N , this condition requires m u < 10 keV Q which is trivially satisfied 
with the observed neutrino mass scale 0.05 eV < m u < 0.33 eV. Similarly, one can find 
that the non-equilibration of the effective operator (2.2) itself is satisfied for a temperature 
below T v : 

— ° i5 -(^T(^r, 

which puts the bound; M < T v , or v s > 1.6 GeV for m u > 0.05 eV and M ~ 10 10 GeV. 

The same argument has to be applied to the supersymmetry breaking operators asso- 
ciated with the supersymmetric operators (2.2). Taking general supersymmetry breaking 
parameters in the scalar potential V, we write 

m u F s 7 _ Km u ~ _ 

V = lhin H Ih^ns, , (3.2) 

v 2 vs v 2 vs 

where Fg is the F-term of S and h. v is a supersymmetry breaking parameter which can 
come from the usual A-term or B-term of the heavy field mass operator, as will be shown 
below. For the moment, we do not specify the origins of these supersymmetry breaking 
parameters keeping their values arbitrary. The requirement of not equilibrating the above 
operators (3.2) leads to the conditions 

— < 5 x 10 7 GeV and — < 5 x 10 5 . (3.3) 

vs v s 
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With these bounds, we find the couplings in Eq. (3.2) satisfying m v Fg/v2Vg < 10~ 5 GeV 
and m v K v /v2Vs < 10~ 7 . These values are too small (compared with the weak scale) to 
lead to any observable consequences on sparticle spectra and couplings probed in collider 
experiments. In particular, it is impossible to realize the mixed sneutrino dark matter 



11, 12] in the framework of Dirac leptogenesis. 



Another important issue in supersymmetric Dirac leptogenesis is the problem with 
unwanted relics. During leptogenesis, the out-of-equilibrium decays of the heavy fields 
($i,$?) produce the scalar components of N and S by the amount of Yx- Since their 
couplings are small as in Eqs. (3.2, 3.3), they never get equilibrated. Thus their initial 
abundances are retained until they decay to the lightest supersymmetric particle (LSP) 
which is considered to be dark matter. If such a non-thermally produced LSP survives, 
the observed dark matter abundance today would put the upper bound 

Y x m x <KT 10 GeV. (3.4) 

This implies Y x < 1(T 12 for the LSP mass around 100 GeV which is clearly in contradiction 
to the condition (2.7) for successful leptogenesis. To avoid this difficulty, we need to have 
large couplings in Eq. (3.2) to make the decay of the scalar n occur while the LSP can 
equilibrate. Requiring the decay temperature To to be larger than the LSP decoupling 
temperature Tlsp, which we take to be 10 GeV, we find 

Similar bounds emerge from the scalar s decay. This consideration is also applicable when 
the scalar n (or s) is the LSP. From Eqs. (3.3, 3.5), the allowed range of the supersymmetry 
breaking parameters is found to be A. v or Fg/vg ~ (10 6 — 10 8 ) GeV. The challenge then is 
to generate such a large supersymmetry breaking parameter. It should be noted again that 
the corresponding contribution to the scalar potential V in Eq. (3.2), being scaled by the 
tiny neutrino mass, is controllably small and is not expected to give rise to any undesirably 
large sparticle mass or any other significant phenomenological consequence at laboratory 
energies. In the next section, we will present a consistent framework for Dirac leptogenesis 
realizing the conditions Fg/vg < TeV but A u ~ 10 GeV with vg ~ TeV. 



4. Origin of vg and A u 

We now show that the conditions for successful Dirac leptogenesis described above can 
be realized consistently in the context of a suitably extended version of nMSSM. As in 
the usual version, let us first couple the singlet superfield S to the Higgs mass operator 
through the superpotential term SH1H2 and thus generate vg ~ TeV dynamically. In this 
scheme, there is no problem with non-thermally produced S as it can equilibrate or decay 
fast by its large coupling with Higgs fields. For a viable scenario of nMSSM containing the 
neutrino sector (2.1), we introduce the following superfields and two global symmetries R 
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and B — L on top of the MSSM: 





S H x H 2 L N $ $ c X 


R 


45 5 18222 
33 3 33333 


B-L 


1 -1 -1 1 


Y 


-f \ -\ 



(4.1) 



Here we have taken the R-charge normalization such that R(W) = 2 and Y as the Standard 
Model hypercharge. This allows us to write down the superpotential 



W 



\SH t H 2 + h ik LH 2 <5> k + ti ik NS$> c k + k x X^ c + ^-X 3 + S . 

3 5 Mp 



(4.2) 



It is now crucial for us to generate high scales for vx = (X) and Fx, which can 
be easily realized by adopting the idea of dynamical supersymmetry breaking fed down by 
gauge messengers || |9|. As an illustration, let us consider the SU(6) x U(l) x U(l) m model 
proposed in Ref. Q. This model contains chiral superfields with the quantum numbers as 
follows: 

.4(15,1,0), JF±(6,-2,±1), 5° (1,3,0), S ±2 (l,3,±2). 

Without contradicting our symmetry assignment (4.1), we can extend the superpotential 
(4.2) to include the terms in Eqs. (2.19), (2.23) and (3.1) of Ref. g. Thus our dynamical 
supersymmetry breaking (DSB) is described by the superpotential 



W, 



DSB 



+ k A X^ 



3 



(4.3) 



where A6 is a Yukawa coupling, Kq is the dynamical scale of SU(6), (j)^ are the messenger 
fields carrying the messenger hypercharge ±1 under U(l) m , and X is the heavy singlet 
superfield used in Eq. (4.2). The strong dynamics of 5^7(6) generates the second (non- 
perturbative) term which drives dynamical supersymmetry breaking in combination of the 
first (tree-level) term. Unlike the original propose for the gauge mediated supersymmetry 
breaking § , we will assume a high scale Aq ~ 10 13 GeV so that the dynamically generated 

1 /2 5/2 1 /2 

F-tevm Fq ~ X^' A 6 ' /Mp' induces rh ~ Fq/Mp ~ TeV, which is a right value for the 
gravitino mass and all the soft masses of the MSSM through gravity mediation. On the 
other hand, dynamical supersymmetry breaking generates a negative mass-squared for the 
J7(l) m -charged scalar field <f> , and thereby vx and Fx are induced through the last two 
terms of Eq. (4.3) as in Ref. |, g. 

Now we can show explicitly how the conditions (3.3, 3.5) are met in our scheme. From 
the superpotential terms (4.2), one finds 



M = k 1 v x , K 



Fx A TP fe3 V X 

— and F s = -r^IW ■ 

vx 5 M% 



(4.4) 



One can easily adjust the Yukawa couplings in Eqs. (4.2, 4.3) to get, for instance, 



v x ~ 10 11 GeV and F x /v x ~ 10 8 GeV , 



(4.5) 
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from which one obtains A u ~ 10 8 GeV and F$ <C m 2 . The latter is too small to become 
relevant in generating the vacuum expectation values of the Higgs bosons and scalar S. 
However, we can have in the scalar potential 



V = h^fs ~ m*s (4.6) 



which play the role of a tadpole contribution in the nMSSM |?]]. This term, together 
with the coupling SH1H2, gives rise to vs ~ TeV and thus appropriate \i and Bfi terms 
respecting the Higgs stability and electroweak symmetry breaking conditions ||. Let us 
finally note that A v is in fact the B-term of the heavy mass operator M <£<I> C . This makes it 
clear why we introduced heavy singlets instead of SU (2) doublets. If the heavy fields ($, $ c ) 
were SU(2) doublets as in Ref. ||, A u would induce via gauge mediation || |[| too large 
soft supersymmetry breaking terms in the MSSM sector characterized by rh ~ a2A ll /4n. 



5. Conclusion 

We have successfully implemented Dirac leptogenesis in an extended version of the nMSSM. 
Unlike that of Murayama and Pierce |§, our extension of the MSSM involves only super- 
fields which are singlets under the Standard Model gauge group. The advantage here is 
the quickening of the otherwise undesirably slow decays of relic singlet sneutrinos, thereby 
avoiding a non-thermal overproduction of LSPs which had plagued earlier Dirac leptoge- 
nesis scenarios. This is achieved by generating heavy mass scales M ~ 10 10 GeV and 
large supersymmetry breaking parameters K u in the range 10 7 — 10 8 GeV by dynamical 
supersymmetry breaking. 

Apart from the usual signatures of the Higgs sector extended with singlet fields, our 
model has no signatures in low energy phenomenology. In particular, the contribution 
of the seemingly large parameter A„ to the scalar potential get scaled by the neutrino 
mass and is inconsequentially small in terms of probing the scalar right-handed neutrino 
sector. A clear laboratory distinction between our model and that for standard Majorana 
leptogenesis is that our light neutrinos are Dirac particles with lepton number conserving 
interactions. Thus, non-observation of neutrinoless nuclear double beta decaying the future 
experiments may hint at the Dirac nature of light neutrinos, providing indirect support for 
Dirac leptogenesis. 
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